Contraints on radiative dark-matter decay from the cosmic microwave background 
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If dark matter decays to electromagnetically-interacting particles, it can inject energy into the 
baryonic gas and thus affect the processes of recombination and reionization. This leaves an imprint 
on the cosmic microwave background (CMB): the large-scale polarization is enhanced, and the small- 
scale temperature fluctuation is damped. We use the WMAP three-year data combined with galaxy 
surveys to constrain radiatively decaying dark matter. Our new limits to the dark-matter decay 
width are about ten times stronger than previous limits. For dark-matter lifetimes that exceed the 
age of the Universe, a limit of (T x < 1.7 x 10~ 2j s" 1 (95% CL) is derived, where £ is the efficiency 
of converting decay energy into ionization energy. Limits for lifetimes short compared with the age 
of the Universe are also derived. We forecast improvements expected from the Planck satellite. 



INTRODUCTION 

After decades of research, the nature of the dark mat- 
ter remains elusive. Nevertheless, many properties of 
dark matter can be inferred from analysis of astrophys- 
ical processes that it might affect. In particular, if the 
dark-matter (DM) particle is unstable, or if is produced 
by decay of an unstable progenitor, then energy released 
during the decay could affect recombination and reion- 
ization. This scenario received much interest [H, 0] af- 
ter the release of the first-year Wilkinson Microwave 
Anisotropy Probe (WMAP) results which showed 
a large temperature-polarization (TE) cross-correlation, 
indicating an earlier epoch of reionization than is eas- 
ily accommodated with the ACDM model In this 
scenario, the decay of the DM particle supplies enough 
energy to provide early reionization. 

Recently, the WMAP team has released data from 
three years of observation and the implications for 
reionization are now more in line with the conventional 
reionization scenario, though with a lowered matter den- 
sity fl m and power-spectrum normalization, the demand 
on early and efficient star formation still exists 0. While 
the original motivation of explaining reionization with 
DM decay is now less compelling, the calculations of the 
effects of DM decay on the CMB can be turned now, with 
the new data, to deriving more stringent constraints to 
decaying DM Q. 

Although 21-cm observations may some day be em- 
ployed to probe particle decay during the dark ages 
the best present constraints come from the CMB. The 
CMB is affected by the ionization history in several ways: 
at large scales, Thomson scattering by free electrons af- 
ter reionization generates polarization from temperature 
anisotropics. This shows up as an enhanced TE cross- 
correlation and polarization auto-correlation (EE). This 
is why the new WMAP data, with greatly improved mea- 
surements of large-angle polarization, will provide much 



more stringent constraints to radiatively decaying DM. 
At smaller scales, Thomson scattering damps the primor- 
dial temperature anisotropy by a factor of e~ 2T '. If the 
matter power spectrum can be measured independently 
by other means, e.g. with large-scale-structure (LSS) sur- 
veys, then this effect can also be used to constrain the 
ionization history. 

In the present work, we use the currently avail- 
able CMB data, including WMAP [1, ACBAR @, 
Boomerang El, CBI [nj], and VSA [12], and the SDSS 
[IH and 2dF [l4j] galaxy surveys, to constrain radiatively- 
decaying DM. This work updates, expands, and improves 
upon our previous work [l| in several ways. In addition 
to using the new WMAP data, which is now far bet- 
ter suited for this analysis than the first-year data, we 
also include LSS data which, as we shall see, provides 
some small additional improvement. We have also im- 
proved our methods, using now Markov Chain Monte 
Carlo (MCMC) techniques for parameter fitting, as we 
have done for analogous constraints to DM annihilation 
in Ref. [HI]. This allows a better understanding of the 
degeneracy between inferred parameters, and it also pro- 
vides more reliable error estimates. This work also differs 
from Ref. [16] by searching a broader range parameter 
space, doing a full likelihood analysis, and including LSS. 
In the regions where our parameter spaces overlap, our 
limit is a factor of a few better than theirs. Refs. [17| and 
[l8t also overlap with this work, but the former studies 
only several particular DM candidates, rather than sur- 
veying the entire decaying-DM parameter space, while 
the latter considers heating/ionization only at recombi- 
nation. 

Our calculation of the effects of DM decay on recom- 
bination, reionization, and heating of the intergalactic 
medium, as well as their subsequent impact on the CMB 
temperature and polarization power spectra follows that 
in Ref. [l] . The ionization rate from DM decay is 
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where T x is the decay width (inverse lifetime) of the 
decaying particle, nib is the mass of the baryon parti- 
cle (hydrogen), (Q x /f2(,) is the initial ratio of DM mass 
to baryonic mass, / is the fraction of the decay energy 
deposited in the baryonic gas, c is the speed of light, 
f x = AE/(m x c 2 ), AE is the energy released in the de- 
cay, Eb = 13.6 eV is the binding energy of hydrogen, and 
t is the cosmological time. The rate at which DM decay 
contributes to the heating of the intergalactic medium 
(IGM) is 



K x = Xhf(fx^x/ n b)mbC 2 T x e' 
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The decay energy deposited in the gas contributes both to 
the ionization and heating of the gas. We model the divi- 
sion between these by an ionizing fraction Xi an d heating 
fraction Xh given by [l[ 



Xi = (1 - x e )/3, Xh = (1 + 2z e )/3, 



(3) 



where x e is the ionization fraction. The effects of decay 
on recombination can be determined completely by two 
independent parameters, T x and £, where ( = // x , and 
is related to the quantity £ defined in Ref. I] by £ = T x (. 

We have modified the recombination code RECFAST 
0, which is used in the CMB code CAMB [H for the 
calculation of the ionization history, to account for these 
extra contributions to ionization and heating. The CAMB 
code is used by the MCMC code COSMOMC as its 
CMB driver. We also modified CDSMDMC so we could vary 
the new parameters T x and £ along with the cosmolog- 
ical parameters. For the present study, we have con- 
sidered the following set of 6 cosmological parameters: 
{flbh 2 , Sl c h 2 , 8, T,n s ,A s }, n s and A s are the the spectral 
index and amplitude of primordial perturbations, and r 
is the optical depth in the absence of DM decay. When 
we include LSS (which also encoded baryon acoustic os- 
cillation information) in our analysis, we also marginalize 
over the galaxy bias b, and use data with k < 0.2Mpc//i 
where the growth is linear. 

When the ionization history is altered, the TE and EE 
spectra are enhanced at large scales (small I). The TT 
spectrum at large I is damped by a factor of e~ 2Ttot , 
where r to t is the total optical depth. If the power- 
spectrum amplitude is not known a priori, the large-^ 
suppression may appear as a low-^ increase. However, 
if the matter power spectrum is obtained from galaxy 
surveys, the degeneracy between rtot and A s can be bro- 
ken. Here we use the LSS power spectrum obtained in 
the SDSS [H and 2dF Q surveys. We will see, though, 
that quantitatively their inclusion improves the decaying- 
DM constraint only a little, as the absolute amplitude is 
not known in this data set, and the constraint comes only 
from the shape of the power spectrum. 

If the decay lifetime is long compared with the age of 
the Universe, the exponential factor in Eqs. fTJ) and @ 
reduce to unity, so we can parameterize the effects of 
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FIG. 1: The marginalized probability distribution function 
(solid curve) and the relative mean likelihood (dashed curve) 
of the (T x parameter in the case of long lifetime, for the 
CMB-only (WMAP 3 year) constraint (green upper curves), 
and the CMB and LSS (SDSS+2dF) constraint (blue lower 
curves). The normalization is such that the maximum of the 
function is 1. 



DM decay fully by the combination CT X - ^ n Fig- [H we 
plot the marginalized probability distribution function 
(PDF) of Cr x and the relative mean likelihood, for both 
the constraint obtained from the CMB only and also that 
obtained by including also LSS. We find that these curves 
peak at (T x = 0, indicating that the CMB is consistent 
with no significant DM decay. The PDF drops rapidly 
at (T x ~ 10 -25 s" 1 , and the 95% limit is at 



Cr x < 2.4 x 10 



-25 _-l 



(4) 



from the CMB only. The constraint is improved slightly 
with the addition of LSS data; the 95% limit is then 



Cr x < 1.7 x 10 



(5) 



in this case. The degeneracy of the bias with the galaxy- 
clustering amplitude limits the power of including LSS. If 
the bias can be fixed by some other measurement, such as 
the three-point correlation function, then there should be 
further improvements better than those we have obtained 
here. For the present data, if we assume C < 0.1, we can 
exclude models with lifetimes T" 1 < 5.9 x 10 23 sec in the 
long-lived case. This improvement over previous result 
l] may provide strong constraint on new decaying-DM 
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models such as those given in Ref. 

For a decay lifetime shorter than the age of the Uni- 
verse at the epoch of reionization (~ 10 16 sec), we need 
to fit the two parameters £ and T x independently. In 
Fig. [21 we plot the excluded region (2a) in the two- 
dimensional parameter space (log 10 C: l°gio -T x ), obtained 
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FIG. 2: The 2cr constraints on the decaying-DM parameter 
space. We show the limits from the WMAP1 analysis of 
Ref. [H, our new constraint from WMAP3+SDSS+2dF, as 
well as our forecast for the reach of Planck. 



after marginalization over the other parameters. The re- 
gion filled with red is excluded by the current CMB and 
LSS data set. The boundary of this region is basically 
a straight line with (T x ^const. The line curves a little 
bit at the short-lifetime end (T x ~ 10~ 13 sec). A fit to 
the boundary curve is given by 




FIG. 3: The 2D contours of the distribution of £P X and the 
background parameters for WMAP3+SDSS+2dF data. The 
color is for the relative mean likelihood 



6.77 + 3.96275a: + 0.25858a: 2 -I- 0.00445a; 3 
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(x > -17), 
(6) 



where a = log 10 T x , and y = log 10 £. We have restricted 
the parameter space to T x < 10 



10 1 

-13 



because for still 



shorter lifetime, much of the decay would happen before 
recombination and would thus not affect the ionization 
history of the Universe. 

We will soon have CMB data of still better quality from 
the Planck satellite [23| , to be launched soon. Anticipat- 
ing this development, we use the Fisher-matrix formal- 
ism [HI, [24| to forecast the future measurement error. 
We assume a sky coverage of 0.65 and adopt a fiducial 
model that best fits the WMAP three- year data. Because 
Planck will determine the B-mode polarization, we eval- 
uate the constraint derived from the TT, TE, EE, and 
BB spectra, and use the covariance matrices of Ref. [lH . 
The resulting constraint is also plotted in Fig. O and a 
fit to it is given by 



3.6 + 2.645a + 0.114a 2 (a > 
-25.531 -a (x<-17), 



-17), 



(7) 



where x = log 10 T x , and y = log 10 (. 

Dark-matter decay may also affect the global fitting of 
cosmological parameters. We expect the decay parameter 
(T x to correlate with optical depth r, as both of these are 
related to the ionization history. We also expect it to cor- 
relate with A s and n s , as scattering from electrons damps 



small-scale CMB anisotropy. There should also be some 
correlation with Q c , as this fixes the rate of input of decay 
energy. We plot the 2D contours of the (T x parameter 
with other cosmological parameters in Fig. [3] There is 
indeed some correlation with the above parameters. By 
contrast, if we consider additional ionization and heat- 
ing due to DM annihilation [l5j (rather than decay) , the 
correlation is very weak because DM annihilation occurs 
mainly during the earlier epoch of recombination when 
the DM density is higher [15]. However, the correlation is 
perhaps not quite as strong as one might have anticipated 
from the Fisher-matrix analysis [l|. This, we believe, is 
because the plotted contours follow marginalization, dur- 
ing which the correlation is weakened somewhat. 

We have studied the impact of decaying DM on the ion- 
ization history and its subsequent effects on the CMB. 
Our investigation is largely model- independent; i.e., it 
does not depend on the particular type of particle, but 
only on the decay width T x and energy-conversion effi- 
ciency £. Dark-matter decay may affect the temperature 
and ionization fraction of the baryonic gas. It may also 
help produce or destroy molecular hydrogen, which plays 
an important role in the formation of the first stars [26j . 
We have not considered such effects in this work, but 
treated DM decay and star formation as two indepen- 
dent processes. We obtained constraints on DM decay by 
using data from WMAP's three-year results and the the 
SDSS and 2dF galaxy surveys. We have derived limits on 
decay parameters [cf., Eqs. (HJ)-©]. The constraint on 



4 



decaying DM will be further improved with future CMB 
experiments such as the Planck satellite [cf., Eq. (JTJ)] , 
and/or direct observation of the dark ages with 21-cm 
surveys [§]. 
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